Two-dimensional spin-filtered chiral network model for the Z 2 quantum spin-Hall 

effect 
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The effects of static disorder on the Z 2 quantum spin-Hall effect for non-interacting electrons 
propagating in two-dimensional space are studied numerically. A two-dimensional time-reversal 
symmetric network model is constructed to account for the effects of static disorder on the propaga- 
tion of non-interacting electrons subjected to spin-orbit couplings. This network model is different 
from past network models belonging to the symplectic symmetry class in that the propagating modes 
along the links of the network can be arranged into an odd number of Kramers doublet, ft is found 
that (f) a two-dimensional metallic phase of finite extent is embedded in a Z 2 insulating phase 
in parameter space and (2) the quantum phase transitions between the metallic and Z 2 insulating 
phases belong to the conventional symplectic universality class in two space dimensions. 

PACS numbers: 73.20.Fz, 71.70.Ej, 73.43.-f, 85.75.-d 



I. INTRODUCTION 

An early triumph of quantum mechanics applied to the 
theory of solids was the understanding that, in the ther- 
modynamic limit, the metallic state can be distinguished 
from the insulating state based on the energy spectrum of 
non-interacting electrons subject to the (static) periodic 
crystalline potential. The Bloch insulating state occurs 
when the chemical potential falls within the energy gap 
between the electronic Bloch bands while the metallic 
state follows otherwise. 

It took another 50 years with the experimental discov- 
ery of the integer quantum Hall effect^ to realize that a 
more refined classification of the Bloch insulating state 
follows from the sensitivity of occupied Bloch states to 
changes in the boundary conditions. A two-dimensional 
electron gas subjected to a strong magnetic field turns 
into a quantum Hall insulating state characterized by a 
quantized Hall conductance in units of e 2 / h££&£&ZJL& 
The topological texture of the quantum Hall insulating 
state manifests itself through the existence of chiral edge 
states^ energy eigenstates that propagate in one direc- 
tion along the boundary of a sample with strip geome- 
try. On the other hand, the topologically trivial Bloch 
insulating state is insensitive to modification of bound- 
ary conditions and, therefore, it does not support gapless 
edge states in a strip geometry. The breaking of time- 
reversal symmetry by the magnetic field in the integer 
quantum Hall effect implies the chirality of edge states: 
all edge states propagate in the same direction. Chiral 
edge states thus cannot be back-scattered into counter 
propagating edge states by impurities. For this reason 
the quantization of the Hall conductance is insensitive 
to the presence of (weak) disorder^ (Strong disorder de- 
stroys the very existence of edge states.) 

The (global) breaking of time- reversal symmetry is not, 
strictly speaking, necessary for integer quantum Hall-like 



physics. As a thought experiment, one can consider, for 
example, a noninteracting two-component electronic gas 
with each component subjected to a magnetic-like field 
of equal magnitude but opposite direction^ Each (inde- 
pendent) component is then characterized by its quan- 
tized Hall conductance. The arithmetic average of the 
two quantized Hall conductances vanishes while their dif- 
ference is quantized in units of 2e* 2 /h with e* the effec- 
tive conserved charge. Bernevig and Zhang in Ref. [ll] 
suggested along these lines that, for some semiconductors 
with time-reversal symmetric noninteracting Hamiltoni- 
ans, the role of the magnetic field is played by the intrin- 
sic spin-orbit coupling while the quantum number that 
distinguishes the two components of the two-dimensional 
noninteracting electronic gas is the electronic spinjii If 
so, the quantized Hall conductance for the electric charge 
(arithmetic average) vanishes while the quantized Hall 
conductance for the spin (difference) is nonvanishing (see 
Fig. [[]). 

In the proposal of Bernevig and Zhang, independent 
quantization of the Hall conductance for each spin re- 
quires two independent U(l) conserved currents. The 
first one follows from charge conservation. The second 
one follows from conservation of the spin quantum num- 
ber perpendicular to the interface in which the electrons 



FIG. I: (Color online:) In the proposal of Ref. [Til , a spin- 
up edge state (full line) and a spin-down edge state (dashed 
line) at the boundary of a two-dimensional electronic droplet 
propagate with opposite velocities. The quantized Hall con- 
ductance for the charge vanishes while that for the spin is 
nonvanishing. 
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FIG. 2: Qualitative plot of the bulk single-particle density 
of state (DOS) as a function of the chemical potential for 
(a) the integer quantum Hall effect without disorder, (b) the 
integer quantum Hall effect with weak disorder, (c) the Z 2 
quantum spin-Hall state without disorder, and (d) the Z 2 
quantum spin-Hall state with weak disorder. The chemical 
potential runs along the vertical axis while the DOS runs 
along the horizontal axis. The gray and white in these fig- 
ures denote extended and localized states, respectively. The 
black straight lines in (b) and (d) denote the critical energies 
at which a quantum phase transition takes place between two 
Hall insulating states for (b) and between a metallic and an 
insulating state for (d). 



are confined. However, while the intrinsic spin-orbit cou- 
pling breaks the SU(2) spin symmetry down to its U(l) 
subgroup, the underlying symmetry responsible for the 
quantization of the spin-Hall conductance in Ref. UlL 
other spin-orbit couplings such as the Rashba spin-orbit 
coupling break this leftover U(l) spin symmetry. This is 
not to say that an unquantized quantum spin Hall effect 
cannot be present if the counterpropagating edge states 
survive the breaking of the residual U(l) spin symmetry. 
However, a physical mechanism different from the one 
protecting the integer quantum Hall effect must then be 
invoked for these edge states to be robust against (weak) 
disorder . 12 - 13 - 14 - 15 - 16 - 17 - 18 - 19 ' 20 - 21 - 22 

Kane and Mele showed in Rcfs. [l?] and [HI that a 
nonintcracting tight-binding Hamiltonian inspired from 
graphenc, with a staggered chemical potential and with 
translation invariant intrinsic and extrinsic (Rashba) 
spin-orbit couplings, realizes a time-reversal symmetric 
insulating state that they dubbed the Z 2 quantum spin- 
Hall state, provided the chemical potential lies within the 
bulk spectral gap (see Fig. Although the SU(2) spin 
symmetry is completely broken in most of coupling space, 
parameter space can nevertheless be divided into two re- 
gions depending on whether the number of Kramers dou- 
blet localized at the edges in a strip geometry is even or 
odd. The dispersion of one Kramers doublet edge state 
must necessarily cross the gap in the bulk of the sample 
when the number of Kramers doublet edge state is odd, 
in which case it supports an intrinsic quantum spin-Hall 
effect: an electric field induces a spin accumulation on 



the edges transverse to the direction of the electric field. 
This insulating state with an odd number of Kramers 
edge state defines the Z 2 quantum spin-Hall state. It 
displays a topological texture different from that of the 
integer quantum Hall statei 18 i 22 ' 23 i 24 The insulating state 
with an even number of Kramers edge state is a conven- 
tional Bloch insulator. 

The effect of disorder is to fill the gap in the bulk 
spectrum of the (clean) Z 2 quantum spin-Hall state. Suf- 
ficiently strong disorder is expected to wash out the Z 2 
quantum spin-Hall state by removing the edge states very 
much in the same way as strong disorder does in the in- 
teger quantum Hall effect. On the other hand, Kane 
and Mele have argued that the Z 2 quantum spin-Hall 
state is robust to a weak time-reversal symmetric disor- 
der, as a single Kramers doublet cannot undergo back- 
scattering by a time-reversal symmetric impurity. Both 
expectations were confirmed by a numerical study of (i) 
the four-probe Landauer-Biittiker and Kubo formula^ 9 - 
and of (ii) the spectral flow induced by changes in the 
twisted boundary conditions ^ By appealing to the sym- 
plectic symmetry of the Kane-Mclc Hamiltonian, tun- 
ing the value of the chemical potential away from the 
tails of the disorder-broaden bands towards their center 
should trigger a disorder-induced transition to a metal- 
lic state (see the two mobility edges below and above 
the metallic state in Fig. [2]).— Onoda, Avishai, and Na- 
gaosa have raised the possibility that the topological na- 
ture of the insulating phase might affect critical prop- 
erties at this transition^ Using standard techniques^ 
to investigate the existence of mobility edges in tight- 
binding Hamiltonian perturbed by on-site disorder (here 
the Kane-Mele Hamiltonian with random on-site ener- 
gies distributed with a box distribution), Onoda et al. 
deduced the existence of a mobility edge separating the 
Z 2 quantum spin-Hall state from a metallic state charac- 
terized by the exponent i/ « 1.6 for the diverging local- 
ization length. This exponent is different from the value 
2-5 < v < 2. 8 28 ' 29 that characterizes the conventional 
mobility edge in the two-dimensional symplectic univer- 
sality class. 

Critical properties at the plateau transition in the in- 
teger quantum Hall effect are the same for two very dif- 
ferent microscopic models. There is the effective tight- 
binding model with random off-diagonal matrix elements 
in the basis of Landau wave functions describing the low- 
est Landau leveh 30 i 31 There is the Chalker-Coddington 
network model valid for disorder potentials that vary 
smoothly on the scale of the cyclotron length.— This 
agreement supports the notion that, for the problem 
of Anderson localization, disorder-induced continuous 
quantum phase transitions fall into universality classes 
determined by dimensionality, intrinsic symmetry, and 
topology. Furthermore, some network models have pro- 
vided useful theoretical insights into the problem of An- 
derson localization and some have even been tractable 
analytically.— The purpose of this paper is to construct 
a network model that realizes a quantum critical point 




FIG. 3: A network model is a collection of equipotential lines 
(of the disorder potential) and nodes. Equipotential lines are 
closed in the bulk and possibly open at the boundaries. In 
this paper, equipotential lines are the boundaries of identical 
squares with rounded corners while nodes are the midpoints 
between adjacent rounded corners. Edge states propagate 
anticlockwise (full lines) or clockwise (dashed lines) along 
equipotential lines if their spin is up or down, respectively. 
Each pair of edge state along an equipotential line can be ar- 
ranged into a single Kramers doublet. Edge states belonging 
two two different equipotential line can exchange their mo- 
menta or spin at the nodes of the network such as S or S'. 
Each node is thus assigned a 4 x 4 unitary scattering matrix. 
Two independent copies of the Chalker-Coddington network 
model are obtained in the limit in which all the 4x4 scatter- 
ing matrices are diagonal with respect to the spin degrees of 
freedom. 



separating the Z 2 quantum spin-Hall state from a metal- 
lic state. From this point of view, the network model 
for the two-dimensional symplectic universality class in- 
troduced in Ref. [H is unsatisfactory as it is built from 
an even number (two) of Kramers doublets propagating 
along the links of the network. Instead, the network 
model that we define in Scc.|n]has a single Kramers dou- 
blet propagating along the links of the network. Spin is 
a good quantum number along the links of the network 
so that the spin-up and spin-down components of the 
Kramers doublet can be assigned opposite velocities (chi- 
ralitics). Scattering takes place at the nodes of the net- 
work. If the scattering matrix is diagonal in spin space, 
the network model realizes the proposal of Bernevig and 
Zhang: two copies of the Chalker-Coddington network 
model for the integer quantum Hall effect arranged so as 
not to break (global) time-reversal symmetry (see Fig. 
[T|). However, we will only demand that the scattering 
matrix at a node respects time- reversal symmetry, i.e., 
it can completely break spin-rotation symmetry. Ran- 
domness is introduced through a spin- independent U(l) 
random phase along the links. We also treat the cases 
of random and non-random scattering matrices at the 
nodes. In either cases, our spin-filtered chiral network 
model captures a continuous quantum phase transition 
between the Z 2 quantum spin-Hall state and the metal- 
lic state. We find in Sec. lIIII the scaling exponent v sa 2.7 
for the localization length that is different from the ex- 



FIG. 4: There are two nonequivalent nodes S and S' from 
the point of view of the transfer matrix. Spin-up (full line) 
and spin down (dashed line) are good quantum numbers on 
the links but need not be conserved by the scattering at the 
nodes. 



ponent v sa 1.6 seen by Onoda et al. but agrees with 
the conventional scaling exponent in the two-dimensional 
symplectic universality class. 

II. DEFINITION 

To represent the effect of static disorder on the co- 
herent propagation of electronic waves constrained to 
a two-dimensional plane and subject to a strong mag- 
netic field perpendicular to it, Chalker and Codding- 
ton introduced a chiral network model in Ref. HH- The 
Chalker-Coddington network model makes three assump- 
tions. The disorder is smooth relative to the character- 
istic microscopic scale: the cyclotron length. Equipoten- 
tial lines of the disorder potential define the boundaries 
of mesoscopic quantum Hall droplets along which chiral 
edge states propagate coherently. Edge states belonging 
to distinct equipotential lines can only undergo a unitary 
scattering process by which momenta is exchanged pro- 
vided the distance between the two equipotential lines is 
of order of the cyclotron length. Such instances are called 
nodes of the network model. 

We are seeking a network model that describes coher- 
ent propagation of electronic waves in a random medium 
that preserves time-reversal symmetry but breaks spin- 
rotation symmetry, in short a symplectic network model. 
A second condition is that the number of edge states that 
propagate along equipotential lines can be arranged into 
an odd number of Kramers doublet. We choose the num- 
ber of Kramers doublet to be one for simplicity. A third 
condition is that the symplectic network model reduces to 
two independent Chalker-Coddington models in some re- 
gion of parameter space. The symplectic network model 
from Ref. [H does not fulfill the last two conditions. 

Given the last condition, it is natural to start with 
spin-filtered edge states moving along equipotential 
lines of the disorder potential depicted as squares with 
rounded corners as is done in Fig. [3J The third condition 
on the symplectic network model is then satisfied when 
all the 4x4 unitary scattering matrices at nodes of the 
network do not couple edge states represented by the ar- 
rows along the full lines with edge states represented by 
the arrows along the dashed lines in Fig. [3] The first two 
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conditions arc otherwise satisfied when all the scattering 
matrices at the nodes of the network model from Fig. [3] 
are the most general 4x4 unitary matrices that respect 
time-reversal symmetry. Without loss of generality we 
choose a node of type S from Fig. [3] The most general 
4x4 unitary scattering matrix that respects time- reversal 
symmetry is given by 

4? 4 ' 

S=( r % tQ ), ° Vf n 

r = tanhx, t = — , 

coshx 

Q = ia cos 9 sin ip l + <r 1 sin 9 cos ip 2 
— a 2 sin sin tp 2 + ct 3 cos 9 cos , 

(2.1) 

with the labeling of incoming and outgoing scattering 
states given in Fig. [U The structure displayed by 
Eq. (|2.1[) can be understood as follows. First, the am- 
plitude for an incoming spin-filtered edge state not to 
tunnel must be spin-independent and thus parametrized 
by the single real number r. Second, the strength of 
quantum tunneling at a node can be parametrized by the 
positive-valued transmission amplitude t that multiplies 
the purely imaginary quaternion Q. The purely imagi- 



nary quaternion acts on the spin- 1/2 degrees of freedom 
through the unit 2x2 matrix a and the Pauli matrices 
a i i 3 an d must therefore depend on four real parame- 
ters. Third, the local gauge transformation S — > W SU 
with U = diag ( e+^i e"^ e +i< ^ e _i ^ ) can absorb 
the dependence of Q on the two independent phase shifts 
< ifi 2 < 27r compatible with time-reversal symmetry. 
At last unitarity delivers the constraints r 2 + 1 2 = 1 and 
QQ> = (T . Up to an overall sign of S and a local gauge 
transformation, S can thus be parametrized by 

{(x,9) | < x < oo, 0<6><7r/2}. (2.2) 

The boundary x = oo for which the transmission am- 
plitude vanishes and the scattering matrix is diagonal 
defines the classical limit of the network model. Quan- 
tum tunneling between neighboring plaqucttes in Fig. [3] 
is very weak when x S> 1. In this limit, the network 
model can be interpreted as follows. The host Z 2 quan- 
tum spin-Hall state, i.e., the translation invariant bulk 
state that supports an odd number of Kramers doublet 
edge states in a confined geometry free of disorder, breaks 
down into droplets of Z 2 quantum spin-Hall states sep- 
arated by smooth random potential barriers. To appre- 
ciate the role played by the parameter < 9 < ir/2, wc 
now consider different values of x and 9 on the boundary 
of parameter space. To this end, it is more convenient to 
replace the scattering matrix by two transfer matrices. 

Nodes of type S in Figs. [3] and 2] are assigned the trans- 
fer matrix M 7 



M = UMU\ 



> -x I s1 

U = diag ( e+iCvi+v^) e -i(^i+v 2 ) e -i(vi-v 2 ) e+iCvi-Va) ) , 

/sinhx coshx sin (9 cos 9 sinhx cos 9 coshx sin 6 

— sin 9 cos 9 sinhx coshx —coshx sin 9 sinhx cos 6 

sinh x cos 9 cosh x sin 9 sinh x cosh x sin 9 cos 9 

\— coshx sin 9 sinhx cos 9 — sin 9 cos 9 sinhx coshx,/ 



M 



cosh 2x — cos 29 



Nodes of type S' in Fig. [3] are assigned the transfer matrix M! 

( \ 

M' 



( ^2 T \ 

V>1 T 



M' = U 1 M'U 2 , 



u, 



M' 



^31 

diag ( e + 5^i +v 2) e-stei+Va) e+Kvi+Va) e -f (vi+y 2 ) ) , 
diag ( e + 2"(^i~^2) e -f (vi-v 2 ) e +l (vi-v 2 ) e -f(vi-^2) ) , 

— cosh x cos sinh x sin sinh x cos 9 — cosh x sin #\ 

— sinhx sin# —coshx cos# coshx sin 9 sinhx cosf 

— sinh x cos 9 cosh x sin cosh x cos — sinh x sin ( 

— cosh x sin — sinh x cos 9 sinh x sin cosh x cos ( 



(2.3a) 



(2.3b) 



FIG. 5: (Color online:) The network model at 9 = de- 
couples into two networks depicted by the thick (red) and 
thin (blue) lines. Full and dashed lines distinguish propaga- 
tion along the links of the networks of up and down spins, 
respectively. 



FIG. 6: (Color online:) The network model at 9 = n/2 de- 
couples into two networks depicted by the thick (red) and 
thin (blue) lines. Full and dashed lines distinguish propaga- 
tion along the links of the networks of up and down spins, 
respectively. 



Here, the convention for initial and final scattering states 
is given in Fig. 0] One verifies that, for all values of < 
ifi 2 < 27r, < x < oo, and < 9 < tt/2 that parametrize 

Ai and M! , the conditions for pseudo-unitary 



A 









-er. 







and time-reversal symmetry 




ia 2 




A 









-icr 2 



= A* 



(2.4) 



(2.5) 



hold (A = M or AC). 

Along the boundary 9 = 0, the transfer matrices (|2.3ap 
and (|2.3b[) reduce to 



M = 



and 



/ cothx 1/ sinh a; 

cothx 1/ sinh a; 

1/ sinh a; cothx 

\ 1/ sinh a; cothx 



(2.6a) 



M' = 




sinh x 

cosh x sinh x 

cosh x 

sinh a; cosh a;/ 



(2.6b) 



respectively. As is depicted in Fig. [5j up and down 
spins have decoupled into two independent Chalkcr- 
Coddington models, each of which describes the integer 
quantum Hall plateau transition. Whenever x — or 
x = oo either M. or M! is diagonal so that edge states 
cannot escape the equipotential lines encircling the local 
extrema of the disorder potential. These are strongly in- 
sulating phases characterized by different integer topo- 
logical (Chern) numbers; one for each spin direction. 
Across the plateau transition the number of edge states 
changes by one for each spin, and so does the number 



of Kramers doublet edge mode. This implies that the 
two insulating phases are distinct in the Z 2 classification. 
Quantum tunneling is strongest at the integer quantum 
Hall transition defined by the condition A4 ~ J\A! for 
which x cc = ln(l + \/2). (By ~ is meant equality in 
magnitude of all matrix elements.) 

Along the boundary 9 = tt/2, the transfer matrices 
(I23a|) and (|2~3b"l) reduce to 



/ 



M 



and 



tanh x 



A I cosh x 





tanh x 
1 / cosh x 






-1/ coshx 
tanhx 




1/ cosha^ 



tanh a; 

(2.7a) 



M' 



sinh x — cosh x\ 

sinh x cosh x 

cosh x — sinh x 

cosh x sinh x 



(2.7b) 



respectively. The network model has decoupled into two 
independent network models as is depicted in Fig. [6] The 
U(l) residual spin- rotation symmetry at 9 = is max- 
imally broken by 9 = tt/2. When x = oo, A4 becomes 
diagonal while M. 1 is off-diagonal so that edge states can- 
not escape the equipotential lines encircling the local ex- 
trema of the disorder potential. The point x = is dom- 
inated by quantum tunneling since M. ~ M 1 arc then 
both anti-diagonal. (By ~ is meant equality in magni- 
tude of all matrix elements.) Furthermore, at x = 0, 
propagation of Kramers doublets is ballistic along de- 
coupled one-dimensional chiral channels. Each network 
model depicted in Fig. [6] belongs to the unitary univer- 
sality class (without topological term) when 9 = tt/2 and 
< x < oo. We thus anticipate an unstable fixed point 
at x = describing a metallic phase and an insulating 
phase for x > 0. 

Along the boundary x = 0, the transfer matrices (|2.3a|) 
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FIG. 7: (Color online:) The network model at x = de- 
couples into two networks depicted by the thick (red) and 
thin (blue) lines. Full and dashed lines distinguish propaga- 
tion along the links of the networks of up and down spins, 
respectively. 



and (|2.3b[) reduce to 



M = 



o cote 

-cote o 

l/sinl 

-l/sin0 



l/sinf 
-l/sin0 

cot ( 

- cot e o 



(2.8a) 



and 



M' 



cos e 








— sin 





— cos e 


sin 6? 








sine 


cos e 





sin 








cos 6 



(2.8b) 



respectively. The network model has decoupled into two 
independent network models as is depicted in Fig. [7] 
When e = 0, Ai is off-diagonal while M. 1 is diagonal 
yielding a strongly insulating phase. The point e = it/2 
is dominated by quantum tunneling since M. ~ M. 1 are 
then both anti-diagonal. (By ~ is meant equality in mag- 
nitude of all matrix elements.) Furthermore, at e = tt/2, 
propagation of Kramers doublets is ballistic along de- 
coupled one-dimensional crural channels. Each network 
model depicted in Fig. [7] belongs to the unitary univer- 
sality class (without topological term) when x = and 
< 9 < tt/2. We thus anticipate an unstable fixed point 
at e = tt/2 describing a metallic phase and an insulating 
phase for < 9 < tt/2. 

Observe that the duality relation 



tanh x = cos 9 



implies that 



M(x = 0,0)~M'(x,0 
M'(x = 0,9) ~ M(x,9 



tt/2), 
tt/2). 



(2.9a) 



(2.9b) 



(By rsj is meant equality in magnitude of all matrix ele- 
ments.) 



From the analysis of the network model on the bound- 
aries of parameter space, we deduce the qualitative phase 
diagram shown in Fig. [8l The numerics of Sec. IIIII con- 
firm the overall topology of this phase diagram. 



The definition of the two-dimensional spin-filtered chi- 
ral network model for the Z 2 quantum spin-Hall effect is 
completed by specifying the boundary conditions. These 
are dictated by the numerical method that we shall use in 
Sec. IIIII Following MacKinnon and Kramer, we seek the 
transfer matrix of a long but narrow sample connected 
at both ends to semi-infinite ideal metallic leads. To 
minimize finite size effects, we impose periodic boundary 
conditions in the transverse direction. The transfer ma- 
trix is then a AM x AM pseudo-unitary matrix that maps 
AM plane waves from the left lead into AM plane waves 
from the right lead that we define as follows. First, we 
consider a slice of the sample that we label by the integer 
n = 1, 2, . . . , N as is depicted in Fig. EJa, (N > AM). We 
assign to this slice the AM x AM pseudo-unitary matrix 
M sl (n), 




FIG. 8: (Color online:) Expected phase diagram from the 
analysis of the network model along the boundaries of param- 
eter space (|2.2|l . The fixed point denoted by a filled (green) 
square along the boundary 6 = is the unstable quantum 
critical point located at x cc = ln(l + v2) separating two in- 
sulating phases in the Chalker-Coddington model. The fixed 
point denoted by the filled (blue) rhombus at the upper left 
corner is the unstable metallic phase. The shape of the metal- 
lic phase is controlled by the symmetry crossover between the 
unitary and symplectic symmetry classes. 
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M sl (n) :=M s ,(n) U™(n) M s (n) vg\n), 

Ujp(ri) =diag( e +^i 1) (") e-^i 1 '^) ■•• e+^W e -#w(«)V 



fM m {n) ••• 
A^i(n) 



^ s (n) 



M 0M (n)\ 






\M M0 (n) 



M M _ x (n) 



(2.10a) 



M MM (n)J 



Ms/(n) = diag(M(n) ■•• KiW)- 

Here, A4 m (n) with m = 1, . . . , M— 1 and A4' m (n) with m = 1, . . . , M are given by Eqs. (|2.3a|) and (|2.3b[) , respectively, 
while we have imposed periodic boundary conditions in the transverse direction with the choice 



M 0Q {n) = M M M(n) 



M 0M (n) = M M o{n) 



2 / sinh x cosh x sin 9 cos 9 

cosh 2x - cos 29 \ - sin 9 cos 9 sinh x cosh x 

2 / sinh x cos sin 9 cosh a; 

cosh 2x - cos 26* \- sin 9 cosh a; sinh x cos 

I 



(2.10b) 
(2.10c) 



The phases </4n (n) with I = 1,2, m = l,...,2Af, and 
n = 1,...,N take values between and 2n. Second, 
we assign to the quasi-one dimensional network model 
depicted in Fig. [SJs the transfer matrix 

N 

M tot := J] M sl (n). (2.10d) 

n=l 

This completes the definition of the two-dimensional 
spin-filtered chiral network model for the Z 2 quantum 
spin-Hall effect. 

We close Sec. [II] by showing that Ai tot belongs to the 
Lie group SO* {AM). By construction, flux conservation, 

•Mtot £ 3 ML = s 3 = fa _° J ® i M , (2.11) 

and time-reversal symmetry, 

Z 2 M tot ^ = M* tot , £ 2 =m 2 ®/ 2M , (2.12) 

where 7 2Af is the 2M x 2M unit matrix hold. It follows 
from Eq. ([2~T2]l that 

Ml t = ^ 2 Ml t ^. (2.13) 
Substituting Eq. (|2~T3|) into Eq. ([2~TT1) yields 

Mot^L^l, 

* = _°>4, (2 ' 14) 

We introduce the matrix 



and write S x = —iAA T . Equation (|2.14p then reads 

^tot ^tat = 1= Mot = AM tot A, (2.16) 

where we have used the identity A 2 = 1. We can rewrite 
Eq. (|2.11[) in terms of M tot , 

M tot E 2 A?Lt = S 2 , 

^ A =(7 1 1 2 )^ i m = ^- (2 - 17) 

With an orthogonal transformation that exchanges rows 
and columns, we can bring S 2 into the form, 

°^° T= {-l M V) =j ™- (2 - i8) 

We thus conclude that O M tot O T is an element of the 
group SO* (AM) defined by the conditions, 

9 J 2 m .9 f = J2M1 9 9 T = I am, 9 G GL(AM, C). 

(2.19) 

III. NUMERICS 

This section is devoted to a numerical study of the 
dependence of the smallest Lyapunov exponent of the 
transfer matrix A4 tot defined in Eq. (|2.10p . as a function 
of the width M of the quasi-one dimensional network 
model. Although Ai tot is taken from a statistical en- 
semble that we will specify below, Lyapunov exponents 
are self-averaging random variables for an infinitely long 
quasi-one dimensional network model, ./V — > ooi 13 i 34 
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n=1 2 ■ ■ ■ N-1 N 



FIG. 9: (Color online:) (a) A slice of the network model assuming periodic boundary conditions in the transverse direction, 
here represented by the gray links, (b) Wire geometry of the network model. 



The eigenvalues of the 4M x AM Hermitian ma- 
trix M. J ot M-tot are doubly degenerate and written as 
exp(±2X,-) with < Xi < X 2 < . . . < X 
calization length £ M is then given by 



The lo- 



M 



lim 



N 



N->oo Xi 



(3.1) 



The localization length £ M is a finite and self-averaging 
length scale that controls the exponential decay of the 
Landaucr conductance for any fixed width M of the in- 
finitely long quasi-one dimensional network model, as the 
transfer matrix (|2.10[) belongs to the group SO*{AM)*2- 
It is of course impossible to study infinitely long quasi- 
one dimensional network models numerically and we shall 
approximate £ M with £ M N obtained from the Lyapunov 
exponents of a finite but iong quasi-one dimensional net- 
work model made of N slices. In our numerics we have 
set N = 5 x 10 5 - 8 x 10 6 . 

As shown by MacKinnon and Kramer j2i criticality in 
two dimensions can be accessed from the dependence of 
the normalized localization length 



A := £ M /M 



(3.2) 



on the width M of the quasi-one dimensional network 
model. For example, if £ denotes the two-dimensional lo- 
calization length and if £ diverges according to the power 
law 



i 



(3.3) 



upon tuning of a single microscopic parameter z close to 
its critical value z c , the singular part of A as M — > oo 
should be given by a scaling function^ 



A~ F(xM 1/v ,CM v ,...). 



(3.4) 



Here, x an d C are the single relevant and dominant irrel- 
evant scaling variables, respectively^ The largest irrele- 
vant scaling exponent satisfies y < 0. We assume that F 
can be expanded in powers of QM V and xM 1 /", 

oo oc 

A ~EE F p/V MW+,/ " (3-5) 

p=0 q=0 

where F p £ R are the expansion coefficients. We also 
assume that the relevant scaling variable \ IS linearly 
related to \z — z c \ while the irrelevant scaling variable £ is 
a constant in the vicinity of the critical point. Finally, for 
any given < 9 < tt/2 from the scattering matrix (|2.1j) . 
wc identify the microscopic parameter z as the parameter 
< x < oo. This motivates the scaling ansatz 



A = E< ) (--4 e) ) 9 ^ 



9 =o 



+ Ea ( S(^-4 9) ) 9 m^ 

q =a 

with the 10 real- valued fitting parameters 



(3.6a) 



and 



v v x {e) A {e) — f {e > 



AS) AS) AS) AS) AS) A6) 

J 0,1' JO, 2; JO, 3' J 1,0' J 1,1' Jl,2- 



(3.6b) 



(3.6c) 



Observe that single-parameter scaling is obeyed by 



oo oo 



A ' : = A - E E F p . g C P X q M p y + ^, (3.7) 

p=l q=0 
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or, in practice, 

A' := A - /{J - 4 9) ) q M V+q,U 

= i2f£ ) q (*-4 9) ) q M^. (3.8) 

9=0 

The values taken by the width M of the quasi-one di- 
mensional network model are M = 4, 8, 16, 32, 64. To 
reduce the statistical error, average over 16 different re- 
alizations of the disorder potential are calculated for any 
given M, x, 6 when 9 is not random and M, x otherwise. 

The disorder potential is modeled in two different ways, 
i.e., we introduce disorder in the transfer matrix (|2.10[) 
as follows. Case I: < < tt/2 and < x < oo are the 
same for all nodes and randomness is introduced by tak- 



1,2, 



in 



1, 



, 2M, 



ing all the phases (jr m (n) with / 
and n = 1, . . . , N to be independently and uniformly dis- 
tributed between and 2tt. Case II: in addition to the 
randomness in the phases 4>m (n) we allow 9 to be in- 
dependently distributed with the probability sin(20) be- 
tween and 7r/2 at each node of the network. As before, 
< x < oo is the same for all nodes. 



A. Case I: Randomness on the links only 

We found two unstable fixed points on the boundaries 
of parameter space (|2 .2[) and the expected phase diagram 
was shown in Fig.[5J Boundary x — oo realizes an insulat- 
ing phase. Boundary 9 = realizes the plateau transition 
at x cc = ln(l + \/2) between two Hall insulating phases in 
the integer quantum Hall effect. Boundaries x = and 
9 = 7r/2 realize the unitary universality class with its in- 
sulating phase terminating at the unstable metallic point 
(x, 9) = (0, ir/2). Any critical point close to the last three 
boundaries is difficult to identify numerically as char- 
acteristic crossover length scales between different uni- 
versality classes become very large. A related difficulty 
comes about from the fact that the characteristic disor- 
der strength can remain stronger than the characteristic 
strength of the spin-rotation symmetry breaking away 
from the boundary 6 = of parameter spaced For this 
reason, we use the scaling ansatz (|3.8[) to search for the 
phase boundaries in the interior of parameter space (|2.2p . 



TABLE I: Critical exponent, normalized localization length, 
and (minimal) node parameter x s as a function of 6. The 
poor agreement at 8 = tt/4 is probably due to the presence of 
a large crossover length scale near the boundary x — 0. 
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V 


A c 




tt/8 
37r/16 
7tt/32 
tt/4 * 


2.85 ±0.30 
2.77 ±0.16 
2.73 ± 0.05 
2.17 ±0.21 


1.87 ±0.09 
1.86 ±0.02 
1.90 ±0.01 
1.82 ±0.01 


0.667 ±0.004 
0.465 ± 0.001 
0.244 ± 0.001 
0.016 ±0.001 




0.0 0.1 0.2 0.3 0.4 

M 1/v lx-x s l 




0.0 0.1 0.2 0.3 0.4 0.5 

M 1/V lx- Xl l 



FIG. 10: (Color online:) (a) Normalized localization length 
A as a function of x at fixed 8 = 37r/16 for the widths M = 4 
(crosses), M — 8 (open squares), M = 16 (filled squares), 
M = 32 (open circles), and 64 (filled circles) of the quasi- 
one dimensional network. Error bars are much smaller than 
symbol sizes. There are two values of x (x s « 0.46 and x l ~ 
0.97) for which A does not appear to depend on M. (b) A fit 
of the data shown in (a) with the help of the one-parameter 
scaling ansatz (|3.8[1 when x is close to i s . Inset: A blow up 
of (a) in the vicinity of x s . (c) A fit of the data shown in (a) 
with the help of the one-parameter scaling ansatz (|3.8|l when 
x is close to x l . Inset: A blow up of (a) in the vicinity of x l . 



For illustration, we present in Fig. [TUTa) the x depen- 
dence of the normalized localization length A for the fixed 
values of 9 = 37r/16 and M = 4, 8, 16, 32. It is seen that 
A increases with increasing M for x between 0.5 and 0.9. 
For fixed 0.5 < x < 0.9, this is cither the signature for 
an extended state or that for a localized state whose lo- 
calization length is larger than the maximal width of the 
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quasi-one dimensional network model. Conversely, for x 
smaller than 0.5 or larger than 0.9, A decreases with in- 
creasing M, i.e., this is the signature of a localized state. 
There appears to be two values of x, that we denote with 
x s < x l , for which A does not depend on M = 4, 8, 16, 32, 
and hence are good candidates for a pair of critical points 
separating a metallic from an insulating phase. The in- 
set of Fig. [T07 b) [Fig. [roT c)] magnifies the dependence 
of A on M = 4,8,16,32,64 close to x s (x t ). On this 
scale x s remains well-defined but not x l . We attribute 
the absence of a single crossing point x l in the inset of 
Fig. \TU[ c) to a large contribution from an irrelevant scal- 
ing variable. This hypothesis is verified in Figs. [TOFb) 
and [TOTc) where the single-parameter dependence of A' 
on the scaling variable M x l v \x — x s \ and M x l v \x — Xj|, 
respectively, is demonstrated (we found the value y ~ — 1 
for the largest irrelevant scaling exponent). The values of 
v, A c , x s , and x l obtained from the scaling ansatz (|3.8p 
for different values of 8 can be found in Tables Q] and [Til 
The values that we obtain for v and A c are consistent 
with those for the standard two-dimensional symplcctic 
universality classi 28 i 29 Our numerical map of the phase 
boundaries separating the metallic from the insulating 
phase in the parameter space (|2.2[) is shown in Fig. [TT] 
The shape of the metallic region in Fig. [11] is controlled 
by the crossover from the unitary to the symplectic sym- 
metry class. 
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0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 
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FIG. 11: (Color online:) Phase diagram for the network 
model in the parameter space (|2.2[1 . The location of a critical 
point denoted by a filled (red) circle follows from the scaling 
ansatz ()3.8|l . That denoted by a cross is a rougher estimate 
due to large symmetry crossover effects. The critical points 
denoted by a filled (green) square and a filled (blue) rhombus 
correspond to the critical points of the Chalker-Coddington 
network model and the unstable metallic fixed point from the 
unitary universality class, respectively. Dashed lines are guide 
to the eye. 



The dependence of the normalized localization length 
A on 8 in the insulating regimes x < ln(l + \/2) and 
x > ln(l + -v/2) are different as is shown in Fig. In the 
small 9 insulating regime x < ln(l±\/2) of Fig.[TlJ A is an 
increasing function of 9 for fixed x and M, as is expected 
from the proximity of a phase boundary to a metallic 
phase. In the insulating regime x > 1 > ln(l + \/2) 
of Fig. [Til A depends weakly on 8 for fixed x and M, 
as is expected from a strongly localized regime. When 
x is held fixed at the Chalker-Coddington critical point 
x cc = ln(l + \/2). A is an increasing function of 9 at 
fixed M while A is an increasing function of M at fixed 
9. This is the expected behavior assuming that any finite 
9 drives the critical point (x, 8) = ( ln(l + V2), 0) into a 
metallic phase. The duality relation (|2.9p is also verified 
numerically in Fig. [T2"l 



B. Case II: Randomness on the links and nodes 



Following Asada et al. in Ref. HjJ we expect that cor- 
rections due to irrelevant scaling variables should be re- 
duced by choosing 8 to be independently distributed be- 
tween and ir/2 with the probability sin(26*) for all nodes 
of the network. As is illustrated with Fig. [13l a metallic 
phase exists when 0.05 < x < 1.0. There are two quan- 
tum critical points x s < x l separating the metallic phase 
x s < x < x l from the insulating phase. The scaling 
analysis must account for an irrelevant scaling variable 
with y w —1 in the vicinity of x s . In the vicinity of x l 
a single-parameter scaling analysis suffices. Both scaling 
analysis, summarized in Table [Till imply that the critical 
points x s w 0.05 and x l w 0.97 belong to the standard 
symplectic universality class. 



TABLE II: Critical exponent, normalized localization length, 
and (maximal) node parameter xi as a function of 9. 



6 


V 


A c 


Xl 


tt/8 


2.78 ±0.19 


1.94 ± 0.10 


0.972 ± 0.004 


37r/16 


2.73 ±0.08 


1.87 ±0.03 


0.970 ± 0.002 


tt/4 


2.65 ±0.06 


1.84 ±0.01 


0.970 ± 0.002 


3tt/8 


2.85 ±0.10 


1.78 ±0.06 


0.982 ± 0.002 



TABLE III: Critical exponent, normalized localization 
length, and node parameter when is distributed between 
and tt/2. 

v Ac Xc 

x s 2.74 ±0.12 1.81 ±0.01 0.047 ±0.001 
x, 2.68 ±0.06 1.82 ±0.01 0.971 ± 0.001 
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FIG. 12: (Color online:) Logarithm of the normalized local- 
ization length A as a function of x at fixed M — 16 for 8 — 
(filled circles), it/ 8 (open circles), 37r/16 (filled squares), n/4 
(open squares), 37r/8 (filled triangles), and tt/2 (open trian- 
gles). The duality relation (|2.9p is verified by plotting the 
dependence of In A on 9 when x — and M = 16 (crosses) 
as a function of x — arctanh(cos 0). Inset: Dependence of A 
on for x = ln(l + \/2) and M = 4 (crosses), M = 8 (open 
squares), M = 16 (filled squares), and M = 32 (open circles). 



IV. SUMMARY 

We have constructed and studied a two-dimensional 
spin-filtered chiral network model for the Z 2 quantum 
spin-Hall effect. Disorder has been implemented in two 
distinct ways. The quantum phase transitions between 
the insulating and metallic state are found to be char- 
acterized by the scaling exponent v ~ 2.7 for the di- 
verging localization length. This value is consistent with 
that found in previous numerical studies of the two- 
dimensional metal to insulator transition in the sym- 
plectic universality class . 28 ' 29 We did not find the value 
v w 1.6, recently observed by Onoda et al. in Ref. |H, 
that was interpreted as the signature of a new univer- 
sality class at the transition between the Z 2 quantum 
spin-Hall insulating and the metallic state. 

It is important to remember the similarities and differ- 
ences between our network model and the lattice model 
studied in Ref. [2^. Common to the two models is that, in 
the absence of disorder, they support a host Z 2 quantum 
spin-Hall state (a host Z 2 insulator for brevity) whereby 
an odd number of Kramers doublet edge states cause 
an accumulation of spin at the edges under an applied 
electric field. The crucial difference lies in the spatial 
correlation of the disorder potential added to the host 
Z 2 insulator. On the one hand, in Ref. [26| disorder is in- 
troduced as a random on-site potential that has no spa- 
tial correlation. On the other hand, our network model 
is obtained by perturbing the host Z 2 insulator with a 
spatially smooth disorder potential that breaks the host 
Z 2 insulator into droplets of Z 2 insulators. The network 



< 



< 



< 




0.2 0.3 

M 1/V lx- Xl l 



FIG. 13: (Color online:) (a) Normalized localization length 
A as a function of x with a random for the widths M = 4 
(times), M = 8 (open squares), M = 16 (filled squares), 
M = 32 (open circles), and 64 (filled circles) of the quasi- 
one dimensional network. Error bars are much smaller than 
symbol sizes. There are two values of x (x g m 0.05 and x t w 
0.97) for which A does not appear to depend on M. (b) A fit 
of the data shown in (a) with the help of the one-parameter 
scaling ansatz (|3.8[) (whereby fq°q — > f q , q for all p,q G N) 
when x is close to x a . Inset: A blow up of (a) in the vicinity 
of x g . (c) A fit of the data shown in (a) with the help of the 



one-parameter scaling ansatz (I3.6a|l (whereby ffi ' — > f Q q and 

fl]i ~ f° r au 9£N) when x is close to x l . Inset: A blow 
up of (a) in the vicinity of x l . 



model can thus be viewed as a coarse-grained effective 
model for Z 2 insulating droplets that are weakly coupled 
through quantum tunneling. 

The intrinsic symmetry (time reversal) respected by 
the statistical ensemble of random Hamiltonians or scat- 
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tcring matrices is not changed by the range of the spatial 
correlation of the disorder. The hypothesis of universality 
would then suggest that the same critical scaling should 
be observed at the localization-delocalization transitions 
in the lattice model of Ref. [26| and in our network model. 
The apparent violation of the universality by the two 
numerical results can be reconciled if one assumes that 
there is a long crossover length scale beyond which micro- 
scopic differences between the two models become irrel- 
evant. Corrections from irrelevant scaling variables may 
strongly depend on the range of the disorder potential, as 
in the case of the plateau transition in the second Lan- 
dau level^ and it could well be that the system sizes 
studied in Ref. [H were not large enough. Verification of 
this scenario is left for future work. 

The fact that our network model is built out of two 
Chalker-Coddington network models coupled in a time- 
reversal invariant way has important consequences. Crit- 
icality in each Chalker-Coddington network model can be 
encoded by the field theory of a single (two-components) 
Dirac fermion coupled to a random vector potential, 
a random mass, and a random scalar potential.— It 



then follows by continuity (for small enough 8) that the 
two lines of critical points emerging from the Chalker- 
Coddington critical point in Fig. [TT] can be encoded by a 
field theory for two flavors of Dirac fermions. It is their 
coupling by disorder that prevents the emergence of a 
time-reversal symmetric and topologically driven quan- 
tum critical bchaviori 40 i 41 
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